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ABSTRACT 

With the advent of very large volume, wide-angle photometric redshift surveys like e.g. 
Pan-STARRS, DES, or PAU, which aim at using the spatial distribution of galaxies 
as a means to constrain the equation of state parameter of dark energy, wde, it has 
become extremely important to understand the influence of redshift inaccuracies on the 
measurement. We have developed a new model for the anisotropic two point large-scale 
(r > 64/i _1 Mpc) correlation function £(r p , 7r), in which nonlinear structure growth and 
nonlinear coherent infall velocities are taken into account, and photometric redshift 
errors can easily be incorporated. In order to test its validity and investigate the 
effects of photometric redshifts, we compare our model with the correlation function 
computed from a suite of 50 large-volume, moderate-resolution numerical iV-body 
simulation boxes, where we can perform the analysis not only in real- and redshift 
space, but also simulate the influence of a gaussian redshift error distribution with an 
absolute rms of a z = 0.015, 0.03, 0.06, and 0.12, respectively. We conclude that for the 
given volume (Vb ox = 2Ah~ 3 Gpc 3 ) and number density (n sa 1.25 x 10~ 4 ) of objects 
the full shape of £(r p ,7r) is modeled accurately enough to use it to derive unbiased 
constraints on the equation of state parameter of dark energy wde and the linear bias 
6, even in the presence of redshift errors of the order of a z = 0.06. 
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1 INTRODUCTION 

T he accelerated expa nsion of the Universe was firs t detected 
bv lRiess et~ai1 (|l99ct ) and lPerlmutter et all j 19991 ) from the 
analysis of supernova type la observations. Over subsequent 
years a variety of independent data sets have confirmed this 
finding dCole et al.l | 2005[ ; IKomatsu et al 112009); Ifficken et al ' 
2009 1; iKessler et aLlhooalSanchez et al.ll2009l ; IPercival eta! 



i.e. a cosmological constant, for which i cqe = —1. (see e.g 



Sanchez etalBoOO; Pcrcival et al. 2010; Bcutler et all 



Blake et al]|201ll; IKomatsu et al1l201ll ; iMehta et all 



20lT|; 



20121 ; 



20101 ; IKomatsu et al.l I2011T ) . The origin of the accelerated 



expansion of the Universe is one the most challenging open 
problems in cosmology. Several theories have been proposed 
to provide an explanatio n for this phenomenon (see e.g. 
iDurrer fe Maartensl 120081 . for a review) . One of the most 
promising solutions is the inclusion of a homogeneously and 
isotropically distributed fluid with a repulsive gravitational 
force to the energy-momentum tensor in the Einstein equa- 
tions. This so-called dark energy can be characterized by 
its equation of state parameter wde = Pde/pde, the ratio 
of the pressure and energy density of this component. Cur- 
rent observations are consistent with the simplest solution, 
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Sanchez et al.ll2012f ). although plenty of room still exists for 
alternative scenarios. 

Several observational probes have been proposed as 
means to obtain new clues about the nature of dark energy. 
Among them, the analysis of Baryonic Acoustic Oscillations 
(BAO s) is considered as on e of the most promising alterna- 
tives (|Albrecht et all 120061 ). BAOs are the fossil signal of 
the acoustic waves that propagated through the photon- 
baryon fluid until the epoch of recombination, when elec- 
trons and protons formed neutral hydrogen and the speed 
of sound dropped rapidly to zero. The maximum distance 
the sound waves were able to travel is called sound horizon. 
This scale is imprinted on the large-scale matter density 
fluctuations. In the two point correlation function £(r) the 
signature of the BAOs s hows up as a sin gle broad bump 
at around 110/i _1 Mpc (jMatsubaral 120041'). in the power 
spect r um as a series of sm all wiggles (jEisenstein fc Hul 
1 19981 ; iMontesano et all 120121 ). Since galaxies form in the 
high-density peaks of the matter density field, BAOs are 
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also present m the distribution ot galaxies at later times , 
and have indee d been detected (lEisenstein et all 2005: 
Cole et al.l 120051; I Cabre fe Gaztan aga 2009]; iBeutler et al.l 



20 111 ; iBlake et alj|201lh locally and at intermediate redshift 



Since they originate in the early Universe where dark energy 
does not play a role and their propagation is described by 
well understood plasma physics, the BAO signal can be used 
as a standard ruler. By measuring the apparent extent of 
the acoustic scale in the directions parallel and perpendic- 
ular to the line-of-sight it is possible to recover the redshift 
evolution of the Hubble parameter H and the angular di- 
ameter distance Da throu gh a simple geometrical relation 
(|Blake fc Glazebrooklliool '). 

Due to the low amplitude of the BAO signal, their anal- 
ysis requires the observation of large volumes, which due 
to the expenditure of time spectroscopy causes, are often 
only accessible with photometric surveys - which in turn 
comes at the price of large redshift inaccuracies. Besides 
this technical difficulty, which we will address in this paper, 
galaxy clustering differs from the linear theory predictions in 
a number of ways that need to be taken into account if these 
measurements are to be used to obtain constraints on cos- 
mological parameters. Nonlinear growth of structure leads 
to coupling between Fourier modes, changing the shape of 
the power spectrum and correlation function. The measured 
clustering statistics are also affected by the gravitationally 
induced peculiar motions of galaxies which introduce a dis- 
tortion when the distance to each galaxy is inferred from its 
observed redshift. Besides this, galaxies are biased tracers of 
the underlying dark matter density field. As a result, the cor- 
relation function of the galaxies could be a modified version 
of that of the mass. Besides their effect on the broad-band 
shape of the correlation function and power spectrum, these 
effects alter the shape and location of t he BAO signature 
(|AnguIo et al.l l2008 : Sanchez et al.ll200Sl ; ISmith et aLlbOOST ) 
and might jeopardize the success of any analysis if they are 
not modelled accurately. 

Redshift-space distortions originating from the pecu- 
liar velocities of the galaxies are unfortunately not the only 
distortions of the line-of-sight component of the correla- 
tion function. The measurement of redshifts are not free 
of errors, which add to the distortion of the clustering sig- 
nal. These errors are particularly important in the case of 
photometric redshifts where they can be as large as cr z = 
0.03 - 0.04, smearing out the clustering signal along the 
line-of-sight, and leading to a significant reduction of its 
amplitude at a given scale. Examples of large photomet- 
ric redshift surveys with the aim to measure the equa- 
tion of state of dark energy using galaxy clustering are the 
Panoramic Su rvey Telescope And Rapid Response Sy stem 
Pan-STARRS (|Chambers fc Pan-STARRS T eam 2004;|), the 



upcoming Dark Energy Survey DES (|Tucker et alj|2010j ). 
or PAU (Physics of the Accelerating Universe, see 
iBenitez et al~l (|2009l )'). 

As described above, in order to obtain constraints on 
both Da(z) and H(z) it is necessary to measure the full 
two-dimensional clustering pattern by splitting the comov- 
ing distance between the galaxies into a component perpen- 
dicular, r p , and parallel, it, to the line-of-sight. Despite this, 
most of the theoretical analyses and observations to date 
have focused on angle-averaged statistics which are sensitive 
instead to the parameter combination D A (z)/H(z). In this 



paper we model the full £(r p ,ir) by taking nonlinear struc- 
ture growth, nonlinear coherent infall velocities, and redshift 
errors into account. The paper is structured as follows: In 
Section[2]the analytic model and the tests performed to com- 
pare it with the results of N-body simulations are discussed 
in detail. The performance of the model for real-, redshift- 
and redshift-error space is described in Section [3] where it 
is also compared to previous work. Finally, in Section [3] the 
most important results of this paper are summarized. 



2 METHOD 

2.1 The L-BASICC II simulations 

The validity of the model we are going to describe in the 
following is tested by comparing it to the anisotropic corre- 
lation function measured in a suite of 50 numerical large- 
volume, medium-resolut ion iV-body simulations, the L- 
BAS ICC II simulations jAngulo et al.l 120081 : ISanchez et all 
l200Sf ) ; we use the snapshots at redshift z = 0.5. The L- 
BASICC II simulations correspond to fifty different real- 
izations of a flat ACDM cosmology with JIm = 0.237, 
n s — 0.954, and erg = 0.77, which is in close agree- 
ment with the latest constraints on cos mological parameters 
from CMB and LSS measurements ( e .g. [ Sanchez et al.ll2006l : 
ISpergel I l2007i : ISanchez et all 12009 ; iKomatsu et all |201ll ). 
Each simulation box covers a volume of (1340 /i _1 Mpc) 3 = 
2.4/i _3 Gpc 3 with a particle mass of Mdm = 1-75 x 
10 12 /i _1 Mo- Haloes were identified by a Friend-of- Friends 
halo finder as gravitationally bound systems with at least 
ten dark matter particles. Hence, the lowest mass halo in 
the catalogue has a mass of M min = 1.75 x 10 13 ft -1 M . A 
de tailed descript i on of the L-BASICC II simulations is given 
in lAngulo et~al] ||200Sf ). 



2.2 The model 

This section describes our model of the anisotropic two-point 
correlation function £(r p ,n) for real, redshift and redshift 
error space. Non-linear stru cture growth, described by t hird- 
order perturbation the ory (|jain fc Bertschingedll994rj . the 
non-linear Kaiser effect {Kaiser 19871; IScoccimarrdl2004l ) and 
redshift errors are taken into account. The inclusion of all 
of these effects results in an accurate description of £(r p , 7r), 
helping to get more reliable estimates of the cosmological 
parameters. 



2.2.1 Non-linear evolution 

The evolved structure of the Universe at lower redshifts is 
not described with sufficient accuracy by linear perturbation 
theory. In f(r) the B AO peak is broadened and shifted to- 
wards smaller scales (lAngulo et al.ll2008l; ISmith et al.l l 2008l : 



C 



rocce fc Scoccimarro I 20081 ; Sanchez et 



alj |2008ft . 



An in- 
correct modelling of this feature will produce biased esti- 
mates of the cosmological parameters. It is particularly im- 
portant to model the nonlinearities correctly if the full shape 
of the correlation function is taken into account in the anal- 
ysis. Nonetheless, non-linearities must also be included even 
if only the BAO peak is used as a standard ruler, as these 
introduce a small shift of its position which can lead to an 
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up to a five times larger shift in the estimated value of wde 
jBlake fc Bridle! 120051 : lEisenstein et al.1 120051 : lAngulo et al.l 
120081 ). Including smaller scales (where nonlinearities are even 
stronger) in the analysis makes an accurate treatment of 
nonlinear clustering growth a necessity. 



Several 
order to 
structure, 



Jjain fc Bertschingerl 



methods have 
describe the 
e.g. third-order 



l200fj. 120091). reno rn 
JCrocce fc Scoccimarrol 




1994 
ized 
l2006al : 



been developed in 
nonlinear growth of 
per turbation theory 
Ijeong fc Komatsul 
perturbation theory 



ICrocce fc Scoccimarro 



iMcDonaldT 120071; Matarrese fc Pietroni 



Bernardeau. Crocce fc Scoccimarrol 

Matarrese fc Pietronil 



perturba tion 



20081 ) and 
theory jMatsubaral 



120081 ; 

La grangi an 



Pietronil |200S|; iTaruva fc Hiramatsul 

Smith. Hernandez- Monteagudo fc Seliakl 



20( 


)8a b 




200S 


2009 



Taruvaetal.1 120091 ; lElia et al.l l20ld ). In our analysis 



we model the non-linear evolution of density fluctuations by 
means of third-order perturbation theory. The main input 
for the evaluation of the third-order perturbation theory 
P(k) is the linear theory powe r spectrum, which we gener- 
ate u sing the CAMB codtQ ( Lewis. Challinor fc Lasenbvl 
120001 ). The numerical routines taken from 



http://gyudon.as.utexas.edu/ komatsu/CRL/index.html 



are then used for the further processing of the linear CAMB 
output. The third-order perturbation theory results in 
terms proportional to the squared linear power spectrum 
which are added to the linear power spectrum; this sum 
is de noted 2 nd order power spectrum (|jain fc Bertschingerl 
1 1994 ). The output obtained by these routines are the 
2 nd order pow er spectra for th e densi ty (Equation (A15) 
and (A16) in iNishimichi et al.1 (|2007l)) and the v elocity 
(Equation (A17) and (A18) in lNishimichi et al.l (|2007l )l field 
as well as the densit y - vel ocity (Equation (A19) and (A20) 
in INishimichi et all (|2007l )) cross power spectrum. These 
are needed for the calc ulation of the non-linear Ka iser effect 
|Scoccimarrdl2004 ). In I Jain fc Bertschingerl (|l994 ) only the 
density - density case is derived, see Equations (17) and 
(18). ' 



2.2.2 From P(k) to £(r P , tt) 

The non-linear correlation function is obtained by Fourier 
transforming the third-order perturbation theory power 
spectrum: £(r) = / P{k)^^k 2 dk. Then, 



assuming 
2 + 7r 2 , we can 



small angles (large redshifts) , such that r 2 
split the distance between galaxies into a component perpen- 
dicular (r p ) and parallel (n) to the line-of-sight, and calcu- 
late the full two-dimensional real-space correlation function 
£re(r P) 7r). 

So far only the clustering of dark matter was consid- 
ered. In real space the clustering of dark matter and col- 
lapsed objects (haloes or galaxies) can be related by the 
bias parameter b, which in the two-point statistics to first 
order approximation can be assumed to be a constant factor 
(|Frv fc Gaztanagalll994 ) at the large scales we are interested 



1 http://camb.info/sources/ 
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Figure 1. Our model of the two-point correlation function of dark 
matter haloes in real space, £ rs (r p , 7r), evaluated at z = 0.5. The 

- 1 Mpc, 



BAO peak is now a ring at about r ■ 



HOfe- 



here corresponding to the dark green feature. The black lines 
indicate the range in which we fit the model: Distances smaller 
than \r\ ^ 64/i -1 Mpc and larger than |r| > 165/i _1 Mpc will 
not be taken into account in the analysis. Only one quadrant is 
used in the fit, since the information in the other three quadrants 
is redundant. The colour bar on the right shows the values of 
£,rs(r p ,n). 



£halo/galaxy{r p ,K) = & 2 £ d m (?> , 7r) 



(1) 



where £,dm{r p ,n) and thaio/gaiaxv(r v ,ir) are the dark mat- 
ter and halo/galaxy ^ r3 (r p ,n), respectively. As can be seen 
in Figure [ll which shows the prediction of the dark matter 
haloes real-space £ rs (r p ,n) = £halo/ galaxy (r P ,ir) evaluated 
at z = 0.5, for the cosmology of the L-BASICC II simula- 
tions, the BAO peak in £(r) is now distributed over a ring 
with a radius of r ~ 110 /i _1 Mpc. In redshift space this sim- 
ple equation does not hold, since the bias parameter also 
changes the size and shape of the redshift space distortions 
(mainly through the quadrupole and hexadecapole contribu- 
tion of the non-linear Kaiser effect, see Equation Q and ([5]) 
below), and is thus not a simple multiplicative factor any- 
more. A precise modelling of the redshift space distortions 
is of vital importance; we will describe how we incorporate 
them in our model in detail in the following section. 



2.2.3 From real to redshift space 

If the exact positions of the objects (galaxies, haloes) were 
known, (, r s(r p ,n) would describe their clustering. However, 
peculiar velocities of these objects lead to a blurring of their 
exact positions. The peculiar velocities are induced by the 
surrounding density field as galaxies move in the local gravi- 
tational potential. There are two different effects, which have 
a different influence on the correlation function of galaxies: 
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on small scales the random motions of galaxies within their 
host dark matter halo make structures appear elongated 
along the line-of-si ght, an effect known as Fingers-of-God 
fFoG |jacksonlll972r i, and lead to a reduction of the cluster- 
ing signal. We ignore the FoGs in this analysis - first of all we 
are interested in very large scales, where their effect is neg- 
ligible, secondly we test our model against the L-BAS1CC 
If halo catalogues, which due to the lack of substructure (or 
galaxies) do not display any FoGs. 

On large scales coherent infall of objects onto large 
structures dominates the redshift - space distortions. This ef- 
fect was first described bv lKaiserl(|l987t ). The formulation is 
based on the distant observer approximation, meaning that 
the structure is distorted in a plane-parallel fashion, and 
therefore only the line-of-sight component is affected. The 
Kaiser effect in its original, linear form i ncludes the linea r 
power spectrum; for details we refer to (|Hamiltonl 1 19921 ). 
In order to obtain unbiased and more accurate constraints 
on th e cosmological parameters, we have to go a step fur- 
ther l|Cole. Fisher fc Weinberg"! 1 1994T). and model the non- 
linear Kaiser effect ( Scoccimarrd rilll) 1 ■ on the anisotropic 
two-point correlation function in redshift space £ zs (r p ,7r). 
Following the theore tical derivation of the linear Ka i ser ef - 
fect for Z za (r p ,ir) in iHamiltor] (|l992T ) and iHamiltonl jl99$ ) 
the anisotropic redshift space correlation function can be 
expressed as 



= Co(r)P (/0 + 6(r)iM/*) + U(r)P*(p), (2) 



where Pi(fi) are the Legendre polynomials and the multi- 
poles £j(r) are given by 



fc(r) = 6 a £«(r) + |/&e«(r) + fftee(r), 



(3) 



1 1 f' 

= ^fb[^e(r)-^e(r)] + ^f[^e(r)-iee(r)], (4) 



(5) 



lfb[ete(r)-£so(r)] + ji 

35- — •2 iBe( - r) -2* 

where £ss{r), £se{r) and fee(r) are the dark matter density- 
density, density-velocity and velocity-velocity correlation 
functions, respectively, Fourier transformed from the cor- 
responding power spectra, and / = is the logarithmic 
derivative of the growth factor D with respect to the scale 
factor of the Universe, which we compute as 



/CO 



n M ■ (i + z) 



(6) 



To allow for the dark energy to be different from a 
cosmological constant (in which case 7 = 0.55), we use 
7 = 0.55 + 0.05 ■ (1 + w DE {z)) for iu D e(z) > -1 and 
7 = 0.55 + 0.02 • (1 + m)de(z)) for wde(z) < -1 (|Linder! 
l2007t ) in all our calculations. 

The barred correlation functions in Equations ([3])— (J5j) 
are given by 



3r~ 



£,se(r')r' 2 dr' 



(7) 
(8) 
(9) 

If the linear power spectrum would be used instead of 



£ee(r) = 3r 3 / £,ee{r')r' 2 dr' 
Jo 

£ee(r) = 5r -5 / £.ee(r')r' 4 dr' 
Jo 
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Figure 2. £ zs (rj,,7r) of dark matter haloes with nonlinear Kaiser 
effect at redshift z = 0.5. As in Figure [T] the black lines indicate 
the range in which we fit the model. 



the 2 nd order power spectra, the derived formulas above 
would result in the linear Kaiser effect. In Figure [2]the con- 
tours of £zs(fp,7r) modeled including the nonlinear Kaiser 
effect are shown. Redshift-space distortions destroy spheri- 
cal symmetry and change the shape and location of the BAO 
ring. The structure, affected by the coherent infall of objects, 
looks squashed along the line-of-sight, which also increases 
the clustering signal compared to £ rs (r p ,7r). In addition to 
the change in the amplitude, redshift-space distortions also 
contain information about the cosmological parameters and 
bias, as both growth rate and linear bias enter into the mul- 
tipoles which describe them (Equations [3HS]). 



2.2.4 From redshift to redshift error space 

So far the model is able to describe redshift space cluster- 
ing as long as the redshifts (which translate into distances 
between the observer and the objects) are determined with- 
out large redshift errors, as in spectroscopic surveys, where 
they are negligible. However, collecting a large number of 
spectroscopic redshifts is time consuming, especially if the 
aim is to observe a large volume at high redshift. A faster 
alternative is to estimate redshifts from deep photometric 
data, which makes it necessary to include the effect of large 
redshift errors on ^ 2S (r p ,7r) in the model. How this can be 
done will be described in the next subsection. 

Traditionally photometric redshifts are derived from ob- 
served fluxes in five or more broad to medium band filters 
l|Baum II1962T ). The probability distribution of the fitted red- 
shifts depends on the spectral type, magnitude and redshift 
of the observed objects, the filter set and the fitting scheme 
(e.g. a neural network or a comparison with a library of 
template spectra, which are used to perform synthetic pho- 
tometry). However, several authors have found the redshift 
error distribution in realistic surveys to be very close to 
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gaussian. For example, ICunha et al. I (|2009l ) calculated the 
full probability distribution function for ~ 78 million SDSS 
DR7 galaxies using photometric observables and weighted 
sampling from a s pectroscopic subsample of the data, and 
ISaglia et al.l (|2012T l estimated photometric redshifts for ob- 
jects found in the Medium-Deep Fields of Pan-STARRSl us- 
ing available spectroscopic surveys (including SDSS spectra) 
as training and/or verification sets. From a direct compar- 
ison between spectroscopic and photometric redshifts ob- 
tained for the same objects, they find the width of the 
distribution of all galaxies in the sample to be er z < 5%, 
with ~ 1% extreme outliers (defined as the fraction of ob- 
jects for which \z pho t - z spec \ > 0.15 x (1 + z spec ) ), and 
<j z /(l + z) = 0.024 and an insignificant fraction of only 0.4% 
outliers for the object types which will eventually be used 
to measure BAOs in Pan-STARRS, namely luminous red 
galaxies (LRGs) at z < 0.5. 

Hence, as a first order approximation we adopt here a 
single gaussian peak to simulate the influence of photomet- 
ric redshift errors on the measurement of £(r p ,7r). Since the 
uncertainty of the redshift estimation mainly results in a dis- 
tortion of the distance along the line-of-sight, redshift errors 
have a similar effect as peculiar velocities, and to first order 
they do not depend on the local density field. 

The main effect of photometric redshift errors is that 
the line-of-sight component of the distance to a galaxy is 
smeared ou t by an amount S x\\. T his se ction follows t he de - 
scription in iBIake fc Brldlel <|2005l ) and IPhleps et all (|2007T ) 
for modeling the redshift errors where the redshift errors 
are assumed to be gaussian distributed such that the spatial 
displacement Sx is given by 



f(8x\\) oc exp 



-0.5 



Sx\ 



(10) 



where a x is the comoving distance corresponding to the 
rms <j z of the redshift error probability distribution func- 
tion at the considered redshift . Equation (|10[) describes the 
probability distribution function for one single galaxy, but 
the correlation function is a pairwise statistic, hence in or- 
der to simulate the impact on the model ^~ er r(r p , tt) we have 
to compute the pairwise error distribution function. For two 
galaxies with errors <5x||,i and 6x\\j the rms of their pairwise 
error would be 



Sxi, 



= J 5x1 ; + 6x?. 



(11) 



If N B galaxies are contained in the survey then there are 



N p = 



N g (N g -l) 



pairs. Every galaxy is contained in N g 



pairs, and the overall convolving function f p (rr) is a sum of 
the Gaussians corresponding to all pairs ij: 



-0.5 



<5.r 



\,ij,n 



(12) 



We estimate this by randomly drawing N g values 6x\\, n from 
Equation (110[) , from which we calculate all possible pairwise 
errors 5x\\ t ij >n . Equation (|12[) is the probability distribution 
of the smearing of the correlation signal along the line-of- 
sight due to the redshift errors. By convolving £ zs (r p ,7r) 
with this function, we model the effect of photometric red- 
shifts on the anisotropic correlation function. The resulting 
£zerr(r p , 7r) for redshift errors with an rms of er z = 0.03 is 
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Figure 3. The model of £zerr(r p ,ir) at z = 0.5 convolved with a 
pairwise redshift error distribution, which assumes gaussian red- 
shift errors with an rms of 0.03. Because the clustering signal is 
now smeared out along the line-of-sight, in order to capture (to 
zeroth order) the same information as in the fit of the correlation 
function in real and redshift space, we replace the spherical shell 
in which we fit the model (see Figures JTJ and J2}) by a cylinder 
(as indicated by the black rectangular box). 



shown in Figure [3] The redshift errors have the same ef- 
fect on large scales as the FoG have on small scales, only 
the order of magnitude is much larger. The clustering sig- 
nal is smeared out and the structure looks extremely elon- 
gated, the BAO ring is almost completely washed out and 
can barely be identified; with increasing size of the redshift 
errors the signal to noise of the BAO ring decreases rapidly, 
until the signal completely vanishes and cannot be used any- 
more to constrain any cosmological parameters. 



2.3 Constraining cosmological parameters 

We will now demonstrate the capability of the model of 
£(r p , tt) described in the previous section to constrain the 
equation of state parameter of dark energy wde and the lin- 
ear bias b, by fitting the correlation function computed from 
the L-BASICC II simulations. The anisot ropic correlation 
functi on is calculated using the estimator of lLandv fc Szalavl 
(|l993l ) in each single box up to a distance of ±300. Oh 1 Mpc, 
and in bins of 3.0/i~ 1 Mpc in both r p and tt, for real, red- 
shift and redshift error space, respectively. In the latter case 
"redshift errors", that is, small offsets, are added to the co- 
ordinate which has been designated the line-of-sight. These 
offsets have been randomly drawn from a gaussian error dis- 
tribution function with a given rms. From the 50 estimates 
of C( 7 "p> 7r ) th e mean is calculated. 

The model (,(r P ,n) i s a l so evaluated from 
(-300/1 -1 Mpc-300/T 1 Mpc) to (300/iT 1 Mpc, 300/T 1 
Mpc), in steps of (0.5/i _1 Mpc,0.5/i _1 Mpc). In order to 
facilitate an accurate comparison with the simulation, 
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which has been calculated with a six times larger bin 
size, we average over the model bins when fitting woe 
and b. While keeping Ho, f2 m and Q,a fixed at the values 
determined by the L-BASICC II simulation, we vary woe 
and b, which both change the amplitude of the correlation 
function and the shape of the redshift space distortions (see 
Section E) . 

When we calculate the correlation function of the sim- 
ulation boxes, we have to assume a fiducial cosmology - for 
simplicity we chose the input parameters of the L-BASICC 
II, i.e. the correct cosmology (and thus the correct cor- 
relation function with the acoustic peak at the position 
predicted by theory). For any other choice of cosmolog- 
ical parameters the measured redshifts and positions on 
the sky translate differently into distances perpendicular 
and parallel to the line-of- sight, which leads to a distortion 
(|Alcock fc Paczvnskil 19791 ). Therefore when we fit the model 
to the data, we have to rescale the distances accordingly: 



and 



D\ (g) 
D{ (z) 



H f (z) 
H* (z) 



(13) 



(14) 



where the quantities for the fiducial cosmology of the model 
are superscripted with / and the "true" cosmology (of the 
simulation) with t, Da(z) is the angular diameter distance 
and H(z) is the evolution factor, 



D A (2 



dz 



l + zj H(z) 



(15) 



find the best fitting values for wbe and b, and to estimate 
their errors. Due to the limited number of realisations it is 
unfortunately not possible to calculate correct (invertible) 
covariance matrices (the full covariance matrix would take 
100 x 100 pixels into account, hence at least 100 x 100 + 1 
independent simulation boxes are required) . If only the vari- 
ance of the correlation function is used when fitting the mean 
£(r p ,7r), the resulting errors of woe and b (in real and red- 
shift space of the order of ~ 2%) are clearly underestimated. 
Instead we assume that the scatter of 50 fits to the single 
realisations can at least partly account for cosmic variance 
and the otherwise ignored contribution of the off-diagonal 
elements in the covariance matrix. Therefore the variance of 
the best-fitting values inferred in this way can be thought 
to represent a more realistic estimate of the errors of the fit 
parameters, while leaving the actual best-fitting values un- 
changed. Indeed, although the resulting values of woe and 
b reassuringly do not depend on the way of fitting, the size 
of the errors does. Throughout this paper we quote errors 
which have been calculated from the variance of the 50 val- 
ues of wbe and b, which are about an order of magnitude 
larger than those inferred from the fit to the mean. 

We restrict the range in which we fit the parameters 
to -1.6 ^ w OE < -0.4 and 1 sC b < 20. 40000 steps are 
sufficient for the chain to converge towards the best fitting 
values and to explore the likelihood. Since the calculation of 
the third order perturbation theory power spectrum is time 
consuming, to speed up the analysis we run the MCMC on a 
grid, for which we construct a library of correlation functions 
f(r) for -1.6 ^ wo e < -0.4 in steps of -0.001. Once the 
random sampling process of the MCMC has chosen a new 
value for woe, the appropriate £(r) for the nearest value of 
woe is selected from the library and used as starting point 
for the calculation of the model (,(r p ,ir) (and woe set to 
w l o%)- 



and 



h (z) = h ■ y (n M (i + zf + n A (i + z ) 3 C!+-d«)) (i6) 

where a is the scale factor (the evolution of a is given 
by the Friedmann equations and today's value is defined 
to be unity) and 6 the angular extent on the sky. The 
model and the halo catalog correlation functions are com- 
pared within a spherical shell of 64ft _1 Mpc ^ ♦/ r| + 7r 2 ^ 
165ft -1 Mpc in both real and redshift space. In redshift 
error space the shape of that shell is distorted due to 
the smearing of the redshift errors, so in order to sample 
similar scales we replaced the spherical shell by a cylin- 
der, as indicated by a rectangle in Figure [3] The cor- 
ners of the rectangle are defined by [64k' 1 Mpc, Oh' 1 Mpc] 
(left lower corner), [64ft _1 Mpc, 165ft -1 Mpc] (left upper 
corner), [l65ft _1 Mpc, Oft -1 Mpc] (right lower corner) and 
[l65ft _1 Mpc, 165ft _1 Mpc] (right upper corner). 

The model was tested for two different cases, in one 
case the information contained in the amplitude was taken 
into account, and in the other case we analytically marginal- 
ized over the amplit ude using the scheme described in 
I Lewis fc Bridle! (|2002l ). In the latter case only the shape of 
the model is examined. 

The fit is performed by means of a Monte Carlo Markov 
Chain (MCMC) after averaging and rescaling in order to 



3 RESULTS 

In this section the results of the fits of £ rs (r p , n) (real space), 
^zs{i~ p ,tt) (redshift space) and £ zerr (r p , it) (redshift error 
space) will be discussed in detail. The analytic model will 
be tested by fitting the equation of state parameter of dark 
energy woe and the linear bias b against the correlation 
function of the halo catalogue from the L-BASICC II simu- 
lations by following the description of Section 12.31 

3.1 Real Space 

In order to investigate the validity of the third order per- 
turbation theory, we first test the model against the dark 
matter correlation function of the L-BASICC II simulations, 
where b = 1.0. The best fitting parameters are found to be 
wo e = -0.992 ± 0.091, b = 0.998 ± 0.073, from which we 
conclude that the nonlinear structure growth is modeled ac- 
curately enough to obtain unbiased estimates of these pa- 
rameters. 

If we now want to describe the clustering of collapsed 
objects like galaxies - or dark matter haloes - which are 
biased tracers of the dark matter density field, we have to 
include the bias (as defined in Equation[l} in our calculation, 
as described in Section [2.2.31 Since in real space the linear 
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bias is only a multiplicative factor which boosts the ampli- 
tude of the dark matter correlation function, but does not al- 
ter its shape, any information about the bias is contained in 
the amplitude at a given radius r = \n% + 7r 2 . Including the 
amplitude in the fit of £ rs (r p , 7r) of the dark matter haloes in 
the 50 L-BASICC II boxes we obtain w OE = -1.010±0.117 
and b — 2.641 ± 0.183, if we disregard the amplitude, we 
cannot fit the bias, but the same value for woe is obtained 
when the amplitude information is disregarded. The mea- 
sured value of wue is in good agreement with the fidu- 
cial value of the simulation. Since the simulation has only 
medium-resolution, it does not contain very small haloes 
(the lowest mass halo has ten particles, which corresponds 
to a minimum mass of M m i n = 1.75 x 10 13 /i _1 Mq), hence 
a relatively large mean bias measured for all dark matter 
haloes is expected. 

In Figure[4]the contours of the best fitting model (black 
solid line) are plotted over the £ rs (r p , tt) of the L-BASICC II 
dark matter haloes the values of £ rs (r p , 7r) are color coded in 
a logarithmic fashion; only the scales taken into account in 
the fit are shown. Figure [5] shows cuts through the £ rs (r p , tt) 
plane along constant r p for both data and best-fitting model 
(including the lcr uncertainty limits calculated from the vari- 
ance of the correlation functions in the single L-BASICC 
II boxes). The deviations between the model and the L- 
BASICC II Simulations are small compared to the errors. 
The obtained value of woe is in good agreement with the 
fiducial value of the simulation. 



3.2 Redshift Space 

In redshift space the exact positions of the galaxies (and 
therefore the correlation function) are distorted due to the 
additional Doppler shift induced by their peculiar veloci- 
ties, and thus mainly affects the line-of-sight components of 
£(r p ,7r). On large (BAO) scales coherent infall dominates, 
which in previous models of the anisotropic correlation func- 
tion was assu med to be lin ear and modeled following the 
desciption in iKaiserl (I1987T ). In this work the non-linear 



shape only 



-1.012 ±0.139 2.518 ±0.646 



shape and amplitude 



-1.020 ±0.147 2.633 ±0.222 



Kaiser effect ( Scoccimarrd 120041 ') is applied t o the model 
of £.zs (r v , tt). The validit y of the approa c h of Scoccimarrol 
(|2004l ) has been tested bv lJennings et al~l (|2011 ) , who found 
a good match to simulations. Since the size and angular de- 
pendence of the effect depends on the bias of the objects (the 
bias is also contained in the quadrupole and hexadecapole 
needed to evaluate the model £ zs (r p , tt), see Section l2.2p . in 
redshift space it is possible to infer the value of b from the 
shape of the correlation function alone, in contrast to real 
space. The results of the MCMC analysis (again fitted in the 
range 64.0 ^ *J r p + tt 2 ^ 165. 0/i -1 Mpc) are summarized in 
Table [T] the comparison of the best fitting model (black 
solid line) with the logarithmically color coded L-BASICC 
II £zs(r p , tt) is shown in Figure [4] 

The correlation function of the L-BASICCS II halo cat- 
alogue can be well described by the model. The model con- 
tours match the L-BASICC II (, zs (r p ,Tr) almost perfectly 
(see Figure [6] and the resulting values of ujde and b are in 
good agreement with the real space estimates, too. As ex- 
pected the errors are larger in redshift space than in real 
space. Also the error on b is larger if the information con- 
tained in the amplitude is ignored, the determined values 



Table 1. w-qe and b from the model £ zs (r p , tt): Mean and vari- 
ance of the 50 L-BASICCS II boxes, fit in the range 64.0 ^ 
\Jr p + tt' 2 ^ 165. Oh -1 Mpc. The fit has been carried out for two 
cases, one where only the shape was input to the fit and the am- 
plitude was marginalized over, and one where both shape and 
amplitude have been taken into account. 



of Woe and 6 are however consistent and do not depend on 
whether the amplitude is taken into account or not. 



3.3 Redshift Error Space 

One of the advantages of using the anisotropic correlation 
function £{r p ,ir) (or its Fourier space equivalent P(fc p ,fc||)) 
to infer cosmological parameters is that in the presence of 
redshift space distortions the clustering measurement per- 
pendicular to the line-of-sight remains almost unaffected, 
while distortions along the line-of-sight can be modeled and 
thus properly taken into account. This makes it a perfect 
tool to use in the case of photometric redshifts, the large 
errors of which lead to a rather dramatic distortion, as ex- 
plained in Section 12.21 In order to investigate the effect of 
photometric redshift errors on the estimate of wue and b, 
we simulated the influence of a gaussian redshift error dis- 
tribution with a rms of a z — 0.015, 0.03, 0.06, and 0.12, 
respectively, on the measurement. 

As described in Section [2. 31 the spherical shell in which 
the fit was carried out has now been replaced by a cylin- 
der (indicated by the rectangular box in the r p ,7r-plane in 
Figure |[3|) in order to compare (to zeroth order) the same 
information. We found that in the last case a large part of 
the clustering signal is smeared out to distances much larger 
than the 300ft _1 Mpc we are calulating our model for, the 
BAO ring disappears, and the noise increases such that an 
accurate estimate of woe and 6 becomes impossible. While 
still fitting £{r p ,n) only up to 165 hT x Mpc, extending the 
model to distances tt = 2000 h~ x Mpc before convolving it 
with the pairwise redshift error distribution allows us to re- 
cover some of the clustering signal for redshift errors (at 
least for a z 0.06), but its information content is limited 
due to the low signal-to-noise of the data on these scales. 
The values of wue and b we find for a z — 0.015, 0.03, 0.06, 
and 0.12 are summarized in Table [2] contours of the corre- 
sponding models of £ zs (r p ,7r) are shown in comparison to 
the logarithmically color coded measurement from the data 
in Figure [7] 

In Figure [HI Figure O Figure [lOl and Figure [11] we show 
again cuts along constant values of r p through £, ze rr(r p , tt) 
for o z = 0.015, 0.03, 0.06, and 0.12, respectively. The model 
fits the data well, although as expected slightly worse than 
in real and redshift space (see Figures [5] and [6| . 

Figure[l2]shows the values of woe and b and their corre- 
sponding errors for the fit of the model to the 50 L-BASICCS 
II boxes including increasing widths of the redshift errors. As 
long as the errors are smaller than er z ~ 0.06, the measure- 
ment is unbiased. The increase of the errors in redshift error 
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r p / [Mpc/h] r p / [Mpc/h] 



Figure 4. Contours from the best fit model plotted over the correlation function calculated from the L-BASICC II dark matter halo 
catalogues in real (left) and redshift space (right), respectively. 
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Figure 5. Cuts through the real space correlation function £rs(r p ,7r) of the L-BASICC II dark matter haloes along fixed r p , black 
solid lines: mean, dotted lines: lcr-deviation calculated from the variance of the 50 boxes, red solid line: best-fitting wCDM model, blue 
dot-dot-dashed line: ACDM case. 
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Figure 6. Cuts through the redshift space correlation function ^2s(^p,7r) of the L-BASICC II dark matter haloes along fixed r p , black 
solid lines: mean, dotted lines: ltr-deviation calculated from the variance of the 50 boxes, red solid line: best-fitting wCDM model, blue 
dot-dot-dashed line: ACDM case. 



WDE 



0.015 


only shape 
shape and amplitude 


-0.965 ±0.298 
-0.980 ± 0.296 


3.660 ± 2.694 
2.704 ±0.333 


0.030 


only shape 
shape and amplitude 


-0.883 ±0.313 
-0.966 ± 0.363 


6.309 ±3.603 
2.622 ±0.415 


0.060 


only shape 
shape and amplitude 


-1.081 ±0.344 
-1.036 ± 0.402 


5.336 ± 3.931 
2.609 ±0.512 


0.120 


only shape 
shape and amplitude 


-1.316 ±0.348 
-1.199 ±0.424 


4.904 ±4.537 
2.295 ±0.554 



Table 2. uidb and b from the fit of the model t;zerr(r p , n) to the 
50 L-BASICCS II boxes for a~ = 0.015, 0.03, 0.06, and 0.12 (from 
top to bottom). Again the fit has been carried out for two cases, 
one where only the shape was input to the fit and the amplitude 
was marginalized over, and one where both shape and amplitude 
have been taken into account. 



space can be expected: Due to the convolution with the pair- 
wise redshift error distribution not only the clustering signal 
is smeared out along the line-of-sight, but so is the noise and 
thus the errors on the cosmological parameters inferred from 
the measurement increase. Therefore, within a given bin in 



the £zerr(r p , 7i~) plane the variance is increased compared to 
redshift space, and the larger the redshift errors, the larger 
the increase. The loss of information contained in the mul- 
tipoles due to the distortion by the convolution with the 
pairwise redshift error distribution function also means that 
the linear bias b cannot be determined accurately by only 
using the shape of £ Z err{r p , n). 

Large redshift errors also increase the probability for 
the MCMC not to converge within the allowed parameter 
space - cosmic variance is still large even in boxes of the 
size of the L-BASICC II Simulations, and in some of them 
the BAO ring is almost invisible even in real space. In such 
cases redshift errors finally destroy all of the information 
that might have been there before, and the fit fails. The 
larger the redshift errors the more catastrophic failures are 
produced. Even if we use the information contained in the 
amplitude, the fraction of boxes where the correlation func- 
tion can not be fit can be as high as ~ 40% for o z — 0.06. 
For a z = 0.12 the model tends to yield biased results, as can 
be seen from the "best-fitting" values. 

The exact size of the redshift error at which the values 
of wde and b can not be measured accurately anymore and 
their errors become unacceptably large certainly depends on 
the exact shape of the redshift error distribution function (a 
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Figure 7. L-BASICC II ^ ze rr{i"p, n) and contours from the model § 2 
(bottom left) and u z = 0.12 (bottom right). 

more Lorentz-like distribution with broad wings will have a 
larger impact on £,{r p ,ir) than a Gaussian with a compara- 
ble width of the core), and it most certainly also depends on 
the volume and/or number density of the survey: since the 
reason the fit fails is mainly that the BAO feature vanishes 
in the increasing noise, the larger the signal-to-noise on large 
scales, the larger the redshift errors may be at which the dis- 
appearance of the BAO ring occurs. There may also be the 
possibility to improve the signal-to-noise ratio for photomet- 
ric data by using the full probability distribution function 
of the redshifts in combination with a set of spectroscopic 
redshifts in the same area and redshift range, a method that 
has been shown to be able to improve the clustering signal 




r p / [Mpc/h] 




-150 -100 -50 50 100 150 

r p / [Mpc/h] 



■r(r p , 7t) for a z = 0.015 (top left), er z = 0.03 ( top right) a z = 0.06 



strength in a man ner equivalent to inc reasing the survey size 
by a factor 4-5 bv lMvers et~aH (|2009l ). 



3.4 Is the projected correlation function w(r p ) an 
alternative? 

For the analysis of the correlation function on small scales 
often the projected correlation function has been calcu- 
lated , which is in t heory independent of any radial distor- 
tions [Peebiel[l980|; [5^isTPeebiei[l98i. For small angles 
r 2 = r p + ty 2 . Thus the projected correlation function is 
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Figure 8. Cuts through the correlation function £ 2 err(?"p, 7r) of the L-BASICC II dark matter haloes along fixed r P7 for redshift errors 
of u z = 0.015, black solid lines: mean, dotted lines: l<T-deviation calculated from the variance of the 50 boxes, red solid line: best-fitting 
wCDM model, blue dot-dot-dashed line: ACDM case. 



defined as 

/oo 
f(r„,7r) cbr. (17) 
- oo 

Note that w(r p ) has dimensions of length. If it were in prac- 
tice possible to integrate to infinity, it would in principle be 
possible to recover the three-dimensional real space correla- 
tion function £(r), and w(r p ) would be far better suited to 
infer cosmological parameters from the spatial distribution 
of galaxies than £(r p ,7r), which suffers from redshift space 
distortions. However, since integrating out even to very large 
distances without significantly increasing the noise is not 
feasible (in particular if the signal is smeared out and the 
amplitude diminished by large redshift errors), the integra- 
tion limits have to b e finite and even rather small, see also 
iNorberg et al. 1120091 for an illustration of this. This means 
that a part of the clustering signal, which depends on the 
pairwise redshift probability distribution function as well 
as on the real and the assumed cosmology, can not be re- 
covered. This is illustrated in Figure 1131 where we show 
w(r p ) for different widths of the assumed redshift errors 
(<r z =0.015, 0.03, 0.06, and 0.12, respectively) and two dif- 
ferent choices of the integration limits, Air = 163.5ft -1 Mpc 
and A-k — 298.5ft" 1 Mpc. The resulting shape of w(r p ) de- 
pends strongly on both the width of redshift errors and the 



size of the integration limits: if the integration limits are 
very large, most of the signal can be recovered and the dif- 
ference between real and redshift space and the correlation 
function affected by errors is, although small, still visible. In 
practice choosing A-k = 298.5ft" 1 Mpc is not advisable, as 
the measurement will be dominated by noise. On the other 
hand, if the BAO ring is supposed to be fully included in 
the integration, the limits can not be much smaller than 
An ~ 150ft" 1 Mpc - in which case the resulting w(r p ) is 
extremely dependent on the size of the redshift errors (i.e. 
the fraction of the signal which can be recovered). 

Figures [141 and [151 show the projected correlation func- 
tion w(r p ) of the L-BASICC II dark matter haloes integrated 
up to TTmax = 298.5ft" 1 Mpc and K max = 163.5ft" 1 Mpc, 
respectively, for redshift errors of a z = 0.015, a z = 0.03, 
a z = 0.06, and a z = 0.12, as well as the best fitting model 
for each case. The amplitude has to be taken into account, 
otherwise the fit fails: there is not enough information in 
the shape alone. The corresponding fitted values of the dark 
energy equation of state parameter wue and the bias b 
and their errors are listed in Table [3] and shown in Fig- 
ure \W\ as a function of a z . Although the fit is not biased for 
a z < 0.06, the errors are, as expected, much larger and more 
quickly increasing with increasing redshift errors than for the 
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Figure 9. As in Figure El but for u, = 0.03. 



»De(tti 



298.5/1- 1 Mpc) 6(tt„ 



298.5/1" 1 Mpc) 3.5 Other work in the literature 



0.0 

0.015 

0.03 

0.06 

0.12 



0.0 

0.015 

0.03 

0.06 

0.12 



-1.018 ±0.326 
-1.018 ±0.419 
-0.980 ± 0.433 
-1.017 ±0.437 
-1.197 ±0.463 



2.477 ±0.717 
2.683 ± 0.551 
2.741 ± 0.551 
2.725 ± 0.608 
2.835 ± 0.862 



WDb(h"ti 



163. 5/1- 1 Mpc) b(ir„ 



163.5k' 1 Mpc) 



-1.078 ±0.386 
-0.962 ±0.363 
-0.979 ±0.372 
-1.034 ±0.420 
-1.143 ±0.576 



2.683 ± 0.470 
2.702 ± 0.446 
2.706 ± 0.456 
2.626 ± 0.550 
2.764 ± 0.804 



Table 3. The best fitting values of wyje and b, as deduced from 
the projected correlation function w(r p ), with integration limits 
of TT m ax = 298.5/i _1 Mpc (top rows) and n ma x = 163. 5/i -1 Mpc 
(bottom rows), respectively, for gaussian redshift errors with cr z = 
0.015, a z = 0.03, a z = 0.06, and a z = 0.12. 



corrsponding fits of 7r). Hence we conclude that for the 
analysis of the large-scale correlation function as a means to 
constrain cosmological parameters from photometric data, 
£(r p ,7r) is better suited than the projected correlation func- 
tion w(r v ). 



In the last decade substantial effort has been invested into 
the accurate modelling of the power spectrum or correla- 
tion function of galaxies, in order to derive tight constraints 
on cosmological parameters. While the available and antici- 
pated data sets have become larger, it has become clear that 
non-linear structure growth, peculiar velocities and galaxy 
biasing have to be described as precisely as possible. A vari- 
ety of different methods have been developed and employed 
to tackle these problems, and the properties of different 
statistics have been investigated. Since ours is the first sys- 
tematic investigation of the influence of redshift errors on the 
measurement of the dar k energy equatio n of state parame- 
ter wue using £(r p ,7r) (|Cai et alj (|2009l 'l have investigated 
the impact of photometric redshift errors on the power spec- 
trum), a direct comparison with other work in the literature 
is not possible. Instead, we will highlight the differences of 
ours to existing models in which redshift errors can poten- 
tially be included. 

iGaztanaga et all (|200ct ) modeled £,(r py tt) in order to es- 
timate the position of the radial acoustic peak = 0h~ 
Mpc, 7r)) and infer the bias b and Om from t he SDSS DR6 
LRG sample (I Adelman- McCarthy et al.ll2008h . For the mod- 
eling (Gaztanagj^^iL 1 20081 ') use linear perturbation theory 
and linear redshift space distortions on large scales, whereas 
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Figure 10. As in Figure [8] but for a = 0.06. 
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our model is non-linear in both the description of structure 
growth and the Kaiser effect. While we do not include Fin- 
gers of God in our model (since they do not occur in a dark 
matter halo catalogue without substructure, such as we use 
to compare our model with), they do consider the one dimen- 
sional velocity dispersion a v , in order to measure /? = O.J r / b 
from it, as discussed in detail in ICabre fc Gaztanagal (|2009l ). 



A second study of (,(r p ,ir) was carried out by 
lOkumura et al.l (|2008l ). Also in their approach both struc- 
ture growth and Kaiser effect were treated line arly, but 



they take the wide angle effect dSzalav et al. 


|l998|) and the 


high-z distortion effect (IMatsubara & Suto! 


1996TI into ac- 


count, which are combined inlMatsubaral (20041). The scales 



examined were set to 60.0h~ < s < 160.0/i~ Mpc, in 
spherical shells like we do in this paper for real and red- 
shift space, similar to our choice. Fitting their model to the 
anis otropic correlation func tion of the SDSS DR5 LRG sam- 
ple (|Eisenstein et al]|2005l ) not only woe and b were con- 
strained but also several other cosmolgical parameters were 
determined. It remains to be tested if taking the wide an- 
gle effect into account improves the accuracy of our model 
such that when the redshift errors become larger than the 
o z > 0.06 at which point our fit becomes extremely inaccu- 
rate, it is still possible to retrieve reliable constraints, but 
this analysis is beyond the scope of this paper). 

In lChuang et~al1 (|2010l ) £(r p , tt) was modeled using halo 



fit (|Smith et alj2003h to take the nonlinear structure growth 
into account. The redshift space distortions are included by 
the linear Kaiser effect on large scales and by t he Fingers 
of God on small scales. In IChuang et ail (|2010l l H(z) and 
Da{z) was constrained instead of wbe- 

Instead of fitting an analytic model to observed data 
it is also possible to fit numerical iV-body simulations (i.e. 
mocks, where the dark matter haloes are populated with 
galaxies either using halo occupation modeling or a semi- 
analytic treatment of galaxy evolution, to which the same 
selection function, mask and survey geometry has been ap- 
plied). Then all kinds of clustering statistics can be calcu- 
lated from the mocks and compared to the ob served data. 
One ex ample for such an attempt is the work of lKazin et al.l 
l|2010al lbh. although they did not measure £(r p , it), but stud- 
ied the detectability of the BAO peak in £(r ) as well as of 
the r adial peak in SDSS LRG DR7 sample (|Zehavi et al. I 
2005). The advantage of using a mock to compare the clus- 
tering statistics with is that nonlinear clustering growth and 
redshift space distortions do not have to be modeled analyt- 
ically, but occur naturally in the simulation. Also if many 
mocks are generated, the covariance matrix can be calcu- 
lated. However, the big disadvantage is of course that a 
proper fit which takes differences in the nonlinear clustering 
growth due to different values of woe into account requires 
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Figure 11. As in Figure[H] but for a z = 0.12. 
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Figure 12. Fitted values of the dark energy equation of state parameter iude (left) and the bias b (right) against the width of the 
redshift errors applied to the L-BASICC II haloes and the model. In black the amplitude and the shape of £(r p , tt) were taken into accout 
for the fit, whereas in red only the shape was considered. 



to cover the full parameter space with a large number of 
mocks, which is extremely time consuming. 

iPadmanabhan &; White! (^OOST ) suggested to utilize the 
multipole moments for estimating cosmological parameters. 
Although this requires a smaller number of mock catalogues 



to calculate the covariance matrix than would be needed for 
£(r p ,7r), it can not be used in the presence of photometric 
redshift errors. 

iKazin et al.l d2Qll j ) transformed the analys i s of th e mul- 
tipole expansion by [Padmanabhan &: White! (|2008 : ) from 
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Figure 14. The projected correlation function w(r p ) of the L-BASICC II dark matter haloes integrated up to 7r max 
for redshift errors of cr z = 0.015 (top right), <r z = 0.03 (top left), <r z = 0.06 (bottom right) and <j z = 0.12 (bottom left), black solid lines: 
mean, error bars: lcr-deviation calculated from the variance of the 50 boxes, red solid line: best-fitting wCDM model, blue dot-dot-dashed 
line: ACDM case, green dotted line: redshift space (<r z = 0.00). 
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Figure 15. The projected correlation function w(r p ) of the L-BASICC II dark matter haloes integrated up to TT max = 163. bhr 1 Mpc 
for redshift errors of <j z = 0.015 (top right), cr z = 0.03 (top left), cr z = 0.06 (bottom right) and tr z = 0.12 (bottom left), black solid lines: 
mean, error bars: lcr-deviation calculated from the variance of the 50 boxes, red solid line: best-fitting wCDM model, blue dot-dot-dashed 
line: ACDM case, green dotted line: redshift space (<r z = 0.00). 
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Figure 16. Fitted values of the dark energy equation of state parameter wy>e (left) and the bias b (right) against the width of the 
redshift errors applied to the L-BASICC II haloes and the model. In black w(r p ) was integrated up to "n ma x = 298. 5/i— IMpc for the fit, 
whereas in red ir max = 163.5/i — 1 Mpc. 
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Fourier space to configuration space. They compared the 
multipole expansion including the hexadecapole with so- 
called "clustering wedges" (£(A M ,s), where /i = sn/s and 
si is the radial component of separation s) to constrain H 
and Da, in order to break the degeneracy between these two 
parameters usually found when only using the monopole. 
The "clustering wedges" are able to provide constraints at 
least at the same level of accuracy as the multipole expan- 
sion. Their fit is based on iV-body simulations, but they also 
argue that for practical use an analytic modeling based on 
physical principles is needed. 

Although our model includes neither the large angle ef- 
fect, nor Fingers of God (which we do not have in the data 
we test it against), since we do take non-linear clustering 
growth and non- linear velocities into account, and we are 
able to fit the linear bias by making use of the informa- 
tion contained in the redshift space distortions, we believe 
it is competitive. Also redshift errors can easily be included, 
which makes it a valuable tool to apply to photometric data, 
if the redshift errors are small enough (a, < 0.06). 



4 CONCLUSIONS 

In this work we developed and tested a model of the 
anisotropic two point correlation function £(r p , n), which we 
used to investigate the influence of photometric redshift er- 
rors on the measurement of the dark energy equation of state 
parameter woe- We modeled £(r p ,7r) using third order per- 
turbation theory |jain fc Bertschingerlll994l ) to account for 
the nonlinear natur e of the growth of s tructure and the non- 
linear Kaiser effect (IScoccimarr ol l2004j V Redshift errors can 
be included in the model by convolving it with the pairwise 
redshift error distribution, which can easily be computed 
from the (known) photometric redshift errors. 

In order to test the validity of our model, we fit it to the 
mean correlation function measured from the dark matter 
haloes in a suite of 50 large-volume, m edium-resolution N- 
body simulations ( the L-BASICCS II, lAngulo et alJl200St 
ISanchez et al.|[2008h . Both in real and redshift space the fit 
yields unbiased values of the dark energy equation of state 
parameter wde and the linear bias b. With approximately 
300 000 haloes per box, in real space wde and b can be deter- 
mined with an accuracy of about 12% and 7%, respectively. 
In redshift space these constraints become slightly weaker, 
woe can be measured with an accuracy of approximately 
15%, and the relative error of b becomes ~ 8%. 

If only the shape information is used to infer woe and 
b, the errors on both will increase due to the lack of infor- 
mation contained in the amplitude. The relative error of the 
bias increases more than the relative error of wde, since the 
value of the bias is mainly encoded in the amplitude (and 
less in the quadrupole and hexadecapole contribution to the 
redshift space distortion), whereas the equation of state pa- 
rameter of dark energy influences both shape and amplitude 
likewise. 

In order to investigate the effect of redshift errors on the 
measurement, we added a small offset to one of the coordi- 
nates of the dark matter haloes, which we drew randomly 
from a gaussian error distribution, and convolved the model 
with the corresponding pairwise redshift error distribution 
in the direction along the line-of-sight (n). Redshift errors 



smear out the clustering signal and diminish its amplitude; 
at the same time the convolution leads to a mixing and in- 
crease of the noise of the measurement in the single pixels, 
because intrinsic errors are also distributed along the line-of- 
sight. The impact of this on the constraints on cosmological 
parameters is two-fold: Since the signal of the BAOs (as the 
main feature of the otherwise smooth correlation function) 
becomes weaker in the observed range of scales, its predic- 
tive power decreases - in the case of very large redshift errors 
(er z > 0.06) the signal is smeared out over such a large range 
of distances that it completely disappears in the noise. How- 
ever, since much higher accuracies can be achieved in realis- 
tic ongoing or nea r-future photomet ric surveys such as e.g. 
Pan-STARRS fsee lSaglia et al.ll2012T l, this is not a cause for 
concern. Int egrating £(r v , n) to obtain w(r p ), as originally 
proposed bv lPeebles (|l980T ) as a means to overcome redshift 
space distortions, does not help to improve the constraints, 
as in real space the BAO is a ring in the £(»>, n) plane, and, 
when integrated, is distributed over 0. < r v < 120/1" 1 Mpc. 
Since it is impossible to integrate £(r p ,n) to n = 00, the re- 
sulting amplitude and shape of w(r p ) depends on the choice 
of integration limits as well as the underlying cosmology, 
which adds a further complication. Secondly, the noise itself 
increases in the presence of redshift errors, which creates an 
additional difficulty. Due to the decreased signal to noise of 
the correlation function, the accuracy of the constraints on 
wde and b decreases. 

In order to beat down systematics coming from cos- 
mic variance (which is still large, even on BAO scales), it 
is desireable (and important) to observe the largest volumes 
possible at one particular redshift. Also, in order to measure 
a possible variation in the equation of state with lookback 
time, observations have to be carried out at higher redshifts 
as well. At this moment in time both is still only feasible with 
photometric redshifts. The anisotropic correlation function 
£(r p , 7r), which can be used to infer cosmological parameters 
like the dark energy equation of state hide, is well suited 
to incorporate photometric redshifts. We have developed a 
model of f (r p , n) which will be able to provide unbiased con- 
straints on wde and b for photometric redshift surveys. The 
maximum redshift error for which this model will work cer- 
tainly depends on the exact shape of the redshfit error dis- 
tribution, the volume and number density of the survey to 
which it is applied. 
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